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ABSTRACT 

SDSS J0924+0219 is the most glaring example of a gravitational lens with anomalous flux ratios: optical 
broad-band photometry shows image D to be a factor of 12 fainter than expected for smooth lens potentials. 
We report spectroscopy showing that the anomaly is present in the broad emission line flux ratios as well. There 
are differences between the emission line and continuum flux ratios: the A/D ratio is 10 in the broad Lyman-a 
line and 19 in the associated continuum. Known variability argues for the presence of microlensing. We show 
that microlensing can account for both the continuum and emission line flux ratios, if the broad emission line 
region is comparable in size to the Einstein radii of the microlenses. Specifically, we need the half-light radius 
of the broad-line region to be Rblr ^ 0.4 Re ~ 9 lt-days, which is small but reasonable. If the broad-line region 
is that large, then stars can contribute only 15-20% of the surface mass density at the positions of the images. 
While we cannot exclude the possibility that millilensing by dark matter substructure is present as well, we 
conclude that microlensing is present and sufficient to explain existing data. Under this hypothesis, the A/D 
flux ratio should return to a value close to unity on a time scale of years rather than millennia. 
Subject headings: cosmology: theory — gravitational lensing — quasars: individual (SDSS J0924+0219) 



1. INTRODUCTION 

The gravitational lens system SDSS J0924+0219 presents 
a fascinating cha llenge for lens modelers. Discovered by 
llnada et al] J2003t) . images from the Sloan Digital Sky Survey 
showed what appeared to be a triple system in a very odd con- 
figuration. The lensing galaxy could be seen in higher resolu- 
tion follow-up observations obtained with the Baade telescope 
at Las Campanas Observatory, which in combination with the 
three known images called for a fourth image. Modeling and 
subtracting the three bright images (A, B, and C) did reveal 
a fourth image, D, but it was a factor of 10 fainter than pre- 
dicted. The image positions are typical of an "inclined quad" 
or "fold" configuration produced by a source near a fold caus- 
tic (see Saha & Williams 2003j), so images A and D ought to 
be nearly equal in brightness (see Keeton et al. 2005). 

Anomalies in lens flux ratios are possible when there 
is small-scale structure in the lensing galaxy, in the 
form of either dark m at ter subhalos (mil li lensin g ; e.g., 
Mao & Schneiderl 1 19981 IMetcalf & MadavJ 1200 ll IChibal 
22Q2|jDalal & Kochanek 200 2]) or stars (microlensing; e.g ., 
Chang & Refsdal 1979; Schechter & Wambsganss! 120021) . 
The challenge is to understand whether either possibility can 
actually explain why image D is so faint. One might hope 
to discriminate between milli- and microlensing by observing 
a component of the QSO that is large compared to the Ein- 
stein rings of stars but small compared t o the Einstein rings of 
subhalos (Moustakas & Metc aill2003t IWisotzki et alJl2003t 
IMetcalf et all2004tlChiba et all2005tlWavth et alJ2005t) . We 
therefore obtained broad emission line flux ratios for SDSS 
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J0924+0219 using the Hubble Space Telescope. We describe 
our observations in §2, our analysis in §3, and our conclu- 
sions in §4. We assume a cosmology with VLm = 0.3, ft\ = 0.7, 
and Ho = 70 k m s" 1 Mpc" 1 . The redshift of the source quasar 
is z s = 1 -524 (jlnada et aD|2003|, and the r edshift of the lens 
galaxy is zi = 0.394 jEigenbrod et al.l2005l) . 

2. OBSERVATIONS 
2.1. Imaging 

SDSS0924 was observed with the Wide Field Channel 
(WFC) of the Advanced Camera for Surveys (ACS) on the 
Hubble Space Telescope, as part of program GO-9744 (PI C. 
Kochanek). Four 547-sec exposures were obtained with the 
F555W filter («V) on 18 November 2003, and four 574-sec 
exposures were obtained with the F814W («I) on 19 Novem- 
ber 2003. An interpolated color composite of the V- and I- 
band images is shown in Figure The four QSO compo- 
nents and the lens galaxy are clearly visible. An Einstein ring 
image of the QSO host galaxy is also apparent in the I-band 
image. 

We performed a non-linear least squares fit to the flattened 
frames from standard ACS data processing, using a photo- 
metric model consisting of four point sources and a de Vau- 
coule urs model galaxy, convolved with a Tiny Tim PSF dKrisd 
119951 v6.2). We fit each exposure separately and used the 
scatter among exposures to assess the measurement uncer- 
tainties. For the galaxy we find an I-band effective radius 
R eS = 0"436 ± Of'004, an axis ratio q = 0.92 ± 0.02, and a po- 
sition angle q = -25° ± 5° (East of North). The positions and 
broad-band fluxes of the QSO components are listed in Ta- 
ble H The HST flux ratios are consistent with those reported 
bv llnada etaTH2003l) . 

2.2. Spectroscopy 

We planned for SDSS0924 to be observed as the final and 
most important target in our HST program GO-9854 to obtain 
spatially resolved spectroscopy of eight quadruply lensed sys- 
tems with the Space Telescope Imaging Spectrograph (STIS). 



FIG. 1 .— Direct ACSAVFC (left) and dispersed ACS/HRC (right) 5" X 5" images of SDSS J0924+0219 from the Hubble Space Telescope. The direct image is 
oriented with North up and East to the left. Clockwise from the top, the quasar components are labeled A, C, B, and D. The dispersed image is rotated by ~70°. 
Each image is displayed with a non-linear stretch given by an inverse hyperbolic sine to highlight both bright and faint features. An approximate wavelength 
scale is superposed on the dispersed image to show the effective wavelength solution. 



Table 1 . Component positions and fluxes of SDSS J0924+02 1 9 
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106.7 ± 1.5 
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-0.4296 


0.0026 
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0.0039 
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NOTE. — Fi and Fv are broad-band fluxes in the HST I and V bands, respectively, expressed in counts per second. The uncertainties derived from the scatter 
among exposures are consistent with Poisson noise. These fluxes can be converted to AB magnitudes using zeropoints of 25.937 in I and 25.718 in V. F\i and F\2 
are spectroscopic continuum fluxes at rest wavelengths of 1700 A and 1216 A, respectively, in units of 10~ 17 ergs cnT 2 s A -1 . Component D has a systematic 
uncertainty of ~20% in F\2 due to galaxy subtraction (see text). F a is the flux in the Lya/N V broad emission line, in units of 10~ 17 ergs cm~ 2 s _1 . 



The observations were attempted on 19 June 2004, but guide 
star acquisition failed. We arranged to reobserve SDSS0924 
with the same setup, but the power supply to STIS was shut 
down before the new observations were carried out. 

We considered whether ACS could be used to conduct the 
desired observations, and realized that the sapphire prism 
(PR200L) in the High Resolution Channel (HRC) of ACS 
would provide coverage of the broad emission lines of 
Lyman-a, Si IV and C IV. The observations were successfully 
executed on 29 May 2005. The combined image of three 880- 
sec exposures is shown in Figure^. The visible and near-UV 
light from 6000 A to 2000 A is dispersed over 80 native HRC 
pixels. The dispersion per pixel ranges from 14000 km s" 1 at 
6000 A to 4000 km s" 1 at 3000 A to 1300 km s" 1 at 2000 A. 

The combined cosmic ray rejected product from STScI pro- 
vided the starting point for spectral extraction. We subtracted 
the galaxy by masking the dispersed footprints of the four 
QSO images and modeling the galaxy su rface brightness dis - 
tribution with a linear b-spline model (cf. Bolton et al. 2005). 
The relati ve wavelength was gi ven by the ground calibration 
of ACS (iPavlovskv et al.1 12004T) . and we shifted the wave- 
length scale of each spectrum to match the obvious emission 
features (Lya/N V, Ly/3/O VI, C IV) present in the QSO dis- 
persed images. The spectrum of each QSO image was con- 



structed by a simple boxcar extraction in each column with a 
boxcar half-width in pixels given by \/A/24. The count rate 
in each spectral pixel was converted to flux using sensitiv- 
ity tables produced from ground calibrations. The extracted 
spectra are shown in Figure|2] 

We fit a simple five parameter model composed of a lin- 
ear continuum and a Gaussian emission line to each spectrum 
over the range 1 120-1460 A in the QSO frame, covering the 
Lya/N V broad emission line. The fits are shown in Figure[2] 
and the line and continuum fluxes are listed in Table [l] as F a 
and F12, respectively. We also report the average continuum 
flux level over the rest wavelength range 1600-1800 A as Fn. 

Figure ^> shows that the spectra for images A and D pass 
fairly close to the center of the galaxy. To assess systematic 
effects due to the galaxy subtraction, we repeated the mea- 
surements with both 5% more and less galaxy subtracted. All 
of the spectral flux measurements in Table ^ varied by less 
than 5%, with the exception of the continuum measurement 
F 12 for image D which varied by 20%. One other possible sys- 
tematic effect is that Fn may contain a small amount (<5%) 
of broad line flux, but that should not affect our conclusions. 

2.3. Comments 
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FIG. 2. — Extracted spectra of the four images, versus observed wave- 
length. Sample errorbars are shown only in certain bins for clarity. The 
vertical lines indicate the expected locations of the broad emission line com- 
plexes Ly ff/O VI, Lyq/N V , C IV, and C III], given the quasar redshift 
7, s = 1.524 ilnada et al. 2003); they provide an external check on our wave- 
length solutions. The red curves show our models of Gaussian emission lines 
at Lya/N V, with linear continua, fit to each spectrum over the rest wave- 
length range 1 120-1460 A. The green lines indicate the average continuum 
levels over the rest wavelength range 1600-1800 A. 



SDSS0924 is anomalous not just in broad-band photometry 
but also in the continuum and broad emission line flux ratios. 
Moreover, the anomalies are different in different passbands: 
A/D is 10 in the emission line, 14-16 in the broad-band fil- 
ters, and >19 in the continuum. That the broad-band values 
lie between the emission line and continuum values makes 
sense if the broad-band filters contain both continuum and 
emission line light. (We estimate that ^10% of the V-band 
flux is from broad emission lines, while the I-band contains 
~20% emission line flux including the 3000 A bump.) For 
comparison, the values for A/C are 5.6 in the emission line, 
2.6-2.9 in the broad-band filters, and >9 in the continuum. 
That the broad-band values lie below both the emission line 
and continuum values may seem puzzling. However, we must 
recall that the photometry and spectroscopy come from dif- 
ferent epochs, an d image C has faded since November 2003 
(Kochanek 2004) such that the broad-band ratio is presently 
A/C~5. 

3. ANALYSIS 
3.1. Basic picture 

Optical broad-band photometry shows that t he images 
vary independently and on a time scale of years (Kochanek 
2004), which is long compared with the predicted time de- 
lays (see §3.2) and therefore implies that microlensing is 
present. Microlensing can also explain why the emission 
line and continuum flux ratios differ, as we shall d emonstrate 
(also seelSchneider & Wam bsganss ll990llMetcalf et alJ2004l 
Wavt h et all2005l) . 

What remains is to explain why the emission line flux ratios 
are anomalous. We initially suspected millilensing, but then 
realized that we would need at least two clumps since both 
C and D have anomalous emission line flux ratios. Since we 
know that microlensing is present, a simpler hypothesis is that 
microlensing produces all of the anomalies. We will show that 
microlensing can indeed explain the data for SDSS0924. 

One obvious question is whether microlensing has the right 
scale to affect broad emission lines. The source plane Einstein 



radius of a star of mass M is 
/ 4GM Di s D os ^ 1/2 
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For comparison, reverberation mapping studies have shown 
that the broad line regions of active galactic nuclei span 
a range of sizes but certainly exte nd down to ^10 lt-days 
and below (e.g., Kaspi et al. 2005, and references therein). 
The comparison is not perfect because reverberation mapping 
studies have focused on Balmer lines (especially H/3) while 
we have observed Lya, and because reverberation mapping 
and microlensing are probably sensitive to source geometry 
in different ways. Nevertheless, it appears that the microlens- 
ing scale is not grossly inappropriate. 

3.2. Macromodel 

For a microlensing analysis we need to estimate the conver- 
gence and shear at the position of each image. We treat the 
lens galaxy as a singular isothermal ellipsoid whose surface 
mass density (in units of the critical density for lensing) can 
be written as 
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where q is the projected axis ratio, q is the orientation angle, 
and b is a mass parameter related to the Einstein radius of the 
lens. We constrain the center of the mass distribution using 
the observed galaxy position. We optimize b and q l but fix 9, 
to ma tch the orientation of the observed galaxy (see Kochane 
2002). We include an external tidal shear to allow for the pos 
sibility that the environment of the lens galaxy (which is cur- 
rently unknown) affects the lens potential (e.g.. lKeeton et alJ 

We fit the image and galaxy positions (not the flux ratios, 
since they are the subject of the microlensing analysis). The 
best-fit model has \ 2 = 33 for v = 2 degrees of freedom, which 
is dominated by the galaxy position suggesting that we may 
have underestimated its uncertainty. (Eliminating the galaxy 
position constraint leads to \ 2 = 2.) Allowing the orientation 
to be free would yield x 2 = 15 for v=\, again dominated by 
the galaxy position, but this model has a misalignment of 48° 
between the mass and the light so we deem it to be implausi- 
ble. 

The model has an axis ratio q mo d = 0.81, which is not very 
different from the observed axis ratio q h s = 0.92. It has an 
external shear 7 ext = 0.10 at position angle 8 ext = 71° (East 
of North); such a shear is typical for four-image lenses and 
sugge sts that the lens may li e in a modest group of galaxies 
(e.g.. lMomcheva et alJl2005l) . The predicted convergence k, 
shear 7, and magnification /i for each image are listed in Ta- 
ble 13 The predicted time delays between the images are all 
less than 12 days. Thus, independent variability in the im- 
ages on a time scale of years cannot be attributed to intrinsic 
variability of the source. 

3.3. Toy microlensing 

A toy model for microlensing features a single star that is 
nearly aligned with one of the lensed images. The image will 
be split into several "microimages" separated by a distance 
comparable to the star's Einstein radius Re- On this tiny scale 
the effects of the galaxy can be approximated as a constant 
convergence n and shear 7, so that the lens equation can be 
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Table 2. Macromodel properties 
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0.560 


-13.4 


D 


0.476 


0.565 


-22.4 



written as 
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in coordinates centered on the star and aligned with the shear. 
For a point source, the magnification of a microimage at (x,y) 
is 



f i(x,y)=Ul- K y-Y 



r e[ r e~ 
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(x 2 +y 2 ) 
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Suppose the star is perfectly aligned with the macroimage, 
so u = v = 0. Consider a negative-parity macroimage, such 
that 1 — k — 7 < and 1 — n +7 > 0. Then the lens equation 
is easily solved to find that there are two microimages at x = 
and y = ±Re / >/l — /e+7, and each has magnification p = 
-1 /[4 7 (1 -K + 7)]. 4 Since the microimages are not separately 
resolved, what matters is t he combined magnification /j, tnt = 
1/[27(1 — K + 7)] (also see ISchechter & Wambsganssl ;2002). 
Using the convergence and shear values from Table |2] yields 
Mtot = -0.813, which corresponds to a flux ratio D/A = 0.031. 

For imperfect alignment, we can still solve the lens equa- 
tion analytically if the source lies on the u-axis or v-axis (see 
Appendix B.l of lKeeton et all2003l for details). For a source 
lying a distance d < R E from the origin along either axis, we 
find D/A = 0.031 + 0.027(<i /R E ) 2 + 0(d /R E ) 4 . 

Clearly microlensing can cause image D to be very faint, if 
there is a star reasonably close to the macroimage. 

3.4. Realistic microlensing 

Going beyond the toy model, we face three important ques- 
tions: (1) When the density of stars is low (so their caustics do 
not intersect), how likely is it that there is a star close enough 
to the macroimage to produce strong flux perturbations? (2) 
As the density of stars increases (so their caustics merge), is it 
still possible to get a strong demagnification? (3) How do the 
conclusions differ for point-like and finite sources? 

To answer these questions, we run microlensing si mulations 
using the ray shooting software by Wambsaanss ( 1990a b, 
1999). We use the convergences and shears from Table |3| 
and assume that a fraction /* of the surface mass density is 
in stars. The software produces a magnification map for each 
image. Treating the continuum region as a point source, we 
can make a histogram of the magnifications in the map to ob- 
tain the probability distribution for the continuum magnifica- 
tion. Given a model of the broad line region, we can convolve 
the magnification map with that model before making the his- 
togram to obtain the probability distribution for the emission 



4 If the macroimage has positive parity (such that 1 - 
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FIG. 3. — Microlensing magnification distributions for the four images of 
SDSS0924, assuming /* = 15% of the surface mass density is in stars. The 
dotted histograms show the distributions for the continuum, assuming it to 
be a point source. The solid histograms show the results for the emission 
lines, assuming the broad line region to be a Gaussian with a half-light radius 
#blr/#£ =0.4. 



line magnification. We treat the BLR as a Gau ssian with half- 
light radius Rblr, since Mortonson et al. (2005) argue that de- 
tails other than the half-light radius do not significantly af- 
fect microlensing magnification distributions. Examples of 
the distributions are shown in Figure [3] 

To analyze both continuum and emission line flux ratios, 
we must determine the joint probability distribution p(/i c ,fii) 
for the continuum and emission line magnifications of each 
image. For each source position, we take /z e from the raw 
magnification map and pi from the convolved map, and then 
use all pixels to construct a histogram in the (p c , ///) plane. 

Finally, to simultaneously analyze images A and D we com- 
pute the four-dimensional probability distribution 5 



p(fi) = p A (n A c ,i4) x p D (p",pf) . 



(6) 



This represents the probability density for image A to have 
continuum and emission line magnifications of [j£ and pf, 
respectively, while image D has /if and pf . The joint four- 
dimensional distribution is just the product of the two two- 
dimensional distributions because microlensing is indepen- 
dent in the two images. 

To quantify microlensing's ability to explain our data for 
SDSS0924, we then compute 



P = J dfip(n)x{ l o 



(4/fJ% > 19 and /if //if > 10 
else 



(7) 



This represents the probability that, if we picked a four-image 
lens at random, it would be at least as anomalous as what we 
have observed. 

This is not really the right figure of merit for evaluating our 
data, though, because we did not select SDSS0924 at random; 
we chose it specifically because it is the most anomalous of 
the 22 known four-image lenses. A better figure of merit is 
the probability of finding at least one strong anomaly in a 



5 To simplify the notation, we collect the four magnifications into the vec- 
tor fi - 
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Fig. 4. — Logarithm of the P- value, or the probability of having a flux ratio 
anomaly more extreme than observed, as a function of the half-light radius of 
the broad line region. Each curve has a fixed value of the stellar mass fraction 
/*; at left, the curves are in order /* = 0.05,0.10,0.15,0.20,0.25 from top 
to bottom. The dotted line shows our nominal acceptability threshold P = 
6.25 X Kr 4 . 



sample of N four-image lenses. Unfortunately, it is not clear 
how to compute P^. For N identical lenses, the binomial dis- 
tribution would give P N = 1 - (1 -P) N w NP (for P < 1). The 
known lenses are not identical, but short of doing microlens- 
ing simulations for all of them (interesting, certainly, but be- 
yond the scope of this work) we must assume that the value 
of P for SDSS0924 is representative. Next, we must deter- 
mine an appropriate value for N. There are 16 four-image 
lenses with "fold" or "cusp" configurations, in which a close 
pair or triplet of images can be analyzed in a robust, model- 
independent way to find flux ratio anomalies ( Keet on et al.l 
I2003L 120051) . We therefore set N = 16, and consider this to 
be somewhat conservative since we are neglecting the non- 
fold/cusp images in these 16 lenses as well as all images in 
the six known "cross" lenses. Lastly, we must decide what 
value of Pn represents a reasonable threshold of acceptability. 
Since we are dealing with complicated posterior probabilities 
we should avoid being too demanding, so we believe it is rea- 
sonable to consider P N = 1%, or P w 0.01 /16 w 6.25 x 10" 4 . 

Figure0]shows P as a function of the half-light radius of the 
broad-line region, for different values of the stellar mass frac- 
tion /» . There are acceptable models with the BLR as large as 
^blr/^£ ~ 0.4, for a stellar mass fraction of /* ~ 15-20%. 
A larger broad-line region is ruled out by these models. A 
smaller broad-line region would allow a wider range of stel- 
lar mass fractions. Thus, we conclude that there is a range 
of microlensing models with reasonable parameters that can 
explain the flux ratio anomaly in SDSS0924 acceptably well. 

There are two limitations to our present analysis. First, we 
assumed for simplicity that the stellar mass fraction is the 
same for images A and D. Second, we assumed that all the 
stars have the same mass. The stellar mass function has very 
littl e effect on the mag nification distribution for a point source 
(see Schechter et al. 2004, and references therein). The same 
result cannot hold for an extended source, because a large 
source must not be sensitive to stars that are sufficiently small, 



but the problem of microlensing of a finite source by unequal- 
mass stars has not yet been studied. Relaxing our two assump- 
tions should only make it easier to find microlensing models 
that fit the data, which means that they do not affect our main 
conclusion that microlensing can provide a reasonable expla- 
nation for SDSS0924. 

4. DISCUSSION 

We have discovered that the continuum and broad emis- 
sion line flux ratios in SDSS0924 differ from each other at 
the factor of two level (for both ima ges C and D). Th at fact, 
together with photometric variability ( Kochanek 2004), estab- 
lishes that microlensing is present in this system. We have 
also found that image D is highly anomalous in both the con- 
tinuum and the broad emission line fluxes. We have shown 
that all of these results can be explained by microlensing. 

The key point is that saddle images can be strongly sup- 
pressed by microlensing. The suppression is generally greater 
for a point source than for an extended source, which is why 
images C and D are fainter (relative to A) in the continuum 
than in the emission line. Even so, microlensing can pro- 
duce a factor of 10 suppression in the emission line flux of 
image D, provided that the QSO broad emission line region 
has a half-light radius Rblr < 0AR E ~ 9 lt-days. While 
this result suggests that /?blr is smaller than we thought 
broad-line regions to be whe n we began our project (also see 
Moustakas & Metcalf 2003), it seems reasonable in light of 
results from reverberat ion mapping of active galactic nuclei 
(e.g.,|Kaspi et al1 l2005l) . That is all we can say without a de- 
tailed understanding of geometric factors that may make sizes 
measured from microlensing different from those measured 
by reverberation mapping. But it is enough for our proof of 
principle that microlensing is sufficient to explain all current 
data for SDSS0924. 

We believe that microlensing offers the best and most nat- 
ural explanation for SDSS0924, but concede that we cannot 
rigorously rule millilensing out. Testing whether millilensing 
is present will require new data, such as flux ratios in narrow 
emission lines (t he narrow line regio n is gen erally th ought 
to be large; e.g.. iKraemer et af]fl9 98: but see B ennert et alJ 
2002 f or a contra sting view) or mid-infrared photometry (e.g., 
Chiba et al. 2005). Another intriguing possibility is spectro- 
scopic variability. If there is only microlensing, then over the 
next few years image D ought to return to the brightness pre- 
dicted by smooth lens models — in both the continuum and 
emission line. If there is any millilensing, then differences be- 
tween observed and smooth model flux ratios will persist for 
centuries. 
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We thank A. Eigenbrod for communicating his measurement 
of the lens galaxy redshift prior to publication. We thank 
Frederic Courbin for valuable comments on the manuscript. 
We thank the referee, Prasenjit Saha, for constructive criti- 
cism and for helpful suggestions about the probability anal- 
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